source("AdvancedFactorFunctions.txt")
data <- read.csv("AthleticsData.csv")
R <- cor(data)
n <- dim(data)[1]
Scree.Plot(R,main="Athletics Data")
FA.Stats(R,n.factors=1:5,n.obs=145, main="RMSEA Plot\n Athletics Data (n=1000)",
         RMSEA.cutoff=0.05)
out <- MLFA(Correlation.Matrix=R,n.factors=3,n.obs=1000,promax.m=3)
Loadings(out,cutoff=.3,num.digits=2)